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We show that graphene-based quantum pumps can tap into evanescent modes, which penetrate 
deeply into the device as a consequence of Klein tunneling. The evanescent modes dominate pumping 
at the Dirac point, and give rise to a universal response under weak driving for short and wide 
pumps, in close analogy to their role for the minimal conductivity in ballistic transport. In contrast, 
evanescent modes contribute negligibly to normal pumps. Our findings add a new incentive for the 
exploration of graphene-based nanoelectronic devices. 



I. INTRODUCTION 

Quantum pumps transfer electrons between two reser- 
voirs by externally varying their scattering properties 
over time. This concept has attracted much attention 
since its inception, due in part to its promise for practical 
applications in nanoelectronics 1 2 and for the definition 
of a current standard, 3 5 but also because of its elegant 
theoretical description in terms of the geometry of the 
control parameter manifold 

Efficient quantum pumping requires strong but energy- 
dependent coupling of the pump to the reservoirs. In nor- 
mal systems, pumping is therefore constrained to propa- 
gating modes, while the poorly coupled evanescent modes 
decay rapidly away from the contacts and therefore can- 
not contribute to the pumped charge. Here, we show that 
the discovery of graphene-based materials^ calls for a re- 
vision of such common concepts in quantum pumping. In 
graphene, the low-energy charge carriers are described by 
a massless Dirac equation, 11 and the unique feature of 
chirality suppresses backscattering at interfaces, result- 
ing in the so-called Klein paradox by which charge carri- 
ers are difficult to confine E2HH1 This seems to inhibit the 
prospects of quantum pumping — unless one properly ac- 
counts for the effects of chirality on the evanescent modes. 
These effects so far have been explored only for station- 
ary ballistic transport, where evanescent chiral electrons 
manifest themselves in macroscopic quantum tunn eling 
close to the Dirac point of charge-neutral graphene! 15 * 16 * 
We show that the scattering of these evanescent modes is 
sufficiently energy-dependent so that they contribute sig- 
nificantly to quantum pumping. Close to the Dirac point, 
they deliver the dominant contribution to the pumped 
charge, which can be characterized by a universal dimen- 
sionless pumping efficiency. 

We start this paper in Sec. [IT] with a description of the 
employed model (a quantum pump driven by two electro- 
static gates) and method (the scattering approach to adi- 
abatic quantum pumping). Our results are presented in 
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Sec. |III| The concluding remarks in Sec. |IV| also contain 
a discussion of practical advantages of graphene-based 
quantum pumps. 
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FIG. 1. (Color online) Left panels: Illustration of a graphene 
quantum pump with highly doped (metallic) contacts. Two 
gates at voltage Vi(t) and V2(t) control the time-dependent 
onsite energy Ui(t) and C/2 (t) in the graphene flake over a 
pumping cycle. This induces a charge transport between two 
contact electrodes, separated by a distance L. The work- 
ing point of the device can be further controlled by an ad- 
ditional back gate with voltage Vq. Right panels: Snapshots 
of the onsite potential along the quantum pump, for highly 
doped (metallic) contacts (top) and undoped contacts (bot- 
tom). The solid base line shows the potential at the working 
point. We show that evanescent modes induced by metallic 
contacts enable finite charge pumping at vanishing nominal 
charge- carrier density, in striking contrast to the case of nor- 
mal pumps where the nanoribbon is replaced by a normal 
conductor. 



II. MODEL AND METHOD 

In order to characterize the unique features of quan- 
tum pumping in graphene, we compare the behavior of 
four different setups: graphene pumps with either heavily 
doped contacts or undoped contacts are contrasted with 
normal-conductor pumps with these two types of con- 
tacts. In the absence of driving (i.e., at the working point 
of the pump), the systems with undoped contacts have 
a uniform carrier density, equal in both the leads and 
in the central pumping region. In contrast, the charge 
carrier density in heavily doped contacts is very large; 
this can be realized by metallic contacts^. We assume 
that the charge carrier density in the system can be fur- 
ther uniformly shifted by a back gate covering both the 
pumping region and leads, thereby tuning the working 
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point of the device. For heavily doped contacts, the back- 
gate-induced change of the charge carrier density in the 
leads can be neglected. Since only the setup with highly 
doped contacts induces evanescent modes, this compari- 
son allows to isolate the requirements for the successful 
deployment of such modes in quantum pumping. 

In all of these four cases, we consider a standard pump- 
ing mechanism, which is shown for the example of a 
graphene pump with heavily doped (metallic) contacts 
in the left panels of Fig. [TJ The pump is operated by 
a cyclic raising and lowering of the potential in two ad- 
ditional electrostatic gates, which modulate the charge 
carrier density in the two halves of the system due to 
their effect on the onsite potential U. This induces charge 
transfer between the two electronic reservoirs of the con- 
tacts. The bottom left panel shows a rotated view of the 
pump, which is further characterized by the width W 
and length L of the pumping region. The top right panel 
shows a snapshot of the onsite potential (dashed line) for 
the case of highly n-doped contacts (the solid line is the 
onsite potential at the working point). The bottom right 
panel shows a corresponding snapshot for undoped leads. 

An elegant formulation of quantum pumping is af- 
forded by the scattering approach^, which considers the 
time dependence of the scattering matrix S(t). In the 
minimal case of adiabatic driving with two independent 
parameters £ = {^1,^2} and single channel reservoirs, 
the charge transferred across the non-interacting scatter- 
ing region reduces to an integral over the area enclosed 
by the driving path in two-dimensional (2D) parameter 
space^, 
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The transmission probability T and the scattering phase 
(j) = a — (3 above are determined by the scattering matrix 
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where r (r') and t (t f ) are reflection and transmission am- 
plitudes for electrons arriving from the left (right) reser- 
voir. 

In a quasi one-dimensional setup with more than 
one channel, indexed by quantum number q, the total 
pumped charge will be a sum over channels Q = ^ Q q 
as long as they can be considered independent. This is 
the case for the quantum pump setup depicted in Fig. [T] 
Across the width W of the system, the onsite potential 
U is y independent, so that different scattering channels 
q remain decoupled. We model this profile by two abrupt 
potential steps in the x direction of equal length L/2 and 
assume the two driving parameters £ = {Ui(i), 
have zero average, and maximum amplitudes 5U\ and 
(5/72, respectively. A back gate is used to control the 
carrier concentration in both pump and leads, which 



we parameterize by the average Fermi momentum kp in 
the pumping region. Undoped contacts have the same 
Fermi momentum kp as the pumping region, while heav- 
ily doped contacts have a much larger Fermi momentum, 
which can effectively be taken as infinite. 19 By employing 
the Dirac equation with negligible inter-valley scattering 
for graphene^, and an effective mass approximation for 
the normal system, this model allows us to compute the 
transmission probabilities and phases T q and <p q by sim- 
ple wave matching. 



III. 



RESULTS AND DISCUSSION 



The finite length L of the pump provides a natural 
scale for the scattering problem. It fixes the energy scales 
E^ = Hvq/L (with Dirac velocity vq) in the graphene 
case, and = h 2 /(2m*L 2 ) (with effective mass m*) 
for the normal conductor, which are related to the level 
spacing of the isolated pump. These energies determine 
two possible driving regimes, depending on the maximum 
amplitude Sui of the dimensionless driving energies ui = 
JJijE^. In the weak driving regime, Sui <C 1, the charge 
pumped in channel q can be approximated by 
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Q q « dZQ q (0) / d Ul du 2 = A u d 2 Q q (0), 



where A u ~ SuiSu2 is the small area enclosed in parame- 
ter space by the driving cycle, wherein d 2 Q q (ui,U2) can 
be approximated by a constant. Away from this regime 
the integral in Eq. ([!]) has to be performed numerically. 

For the graphene pump with heavily doped contacts, 
Eq. ([2} yields 
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sin 2 (k x L) [sm(2k x L) — 2k x L cos(2k x L)] 
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Here the ± sign denotes whether the pump i s doped 
with electrons (plus) or holes, and k x = \/k 2 F — q 2 
is the electron's momentum along the transport direc- 
tion in the pumping region. Momentum k x is real for 
propagating modes |g| < kp and imaginary for modes 
\q\ > kp that are evanescent in the pump. In contrast, 
a graphene pump with undoped leads has no incoming 
lead modes that become evanescent in the pump (in the 
limit Sui C 1), since the Fermi momentum in the pump 
at the working point is identical to the Fermi momentum 
in the contacts. A graphene setup with undoped leads 
pumps a charge 
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kpL 2(qL) 2 cos 2 (k x L) sin 3 (k x L) 
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where k x is constrained to real values since |g| < kp- 

In both cases the pumped charge has a prefactor A u kp, 
indicating that pumping is proportional to the pump's 
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FIG. 2. (Color online) Momentum distribution of pumped 
charge per mode \q as a function of mode index q for varying 
carrier concentration (parameterized by the Fermi momentum 
Blue and brown represent opposite directions of pump- 
ing (left to right or right to left). In graphene, the propagating 
mode with q = (normal incidence) cannot be pumped due 
to the Klein paradox. In the case of graphene with heavily 
doped leads, however, significant pumping is possible due to 
the contribution of the evanescent modes (\q\ > /cf, delineated 
by the dashed line), which dominate around the Dirac point 
(kF — 0). The other pumps can only drive current through 
the propagating modes. 



number N p = AkrW/n of propagating electrons at the 
Fermi energy and the dimensionless driving strength A u . 
By factoring out these two quantities, we obtain the di- 
mensionless pumping response 



d 2 
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which depends only on the system's scattering charac- 
teristics at a given energy. Summing over all incoming 
modes gives the fraction of (propagating) electrons in the 
pump that are pumped per cycle and per unit dimension- 
less driving strength A u , 



(6) 



This defines the total dimensionless pumping response 
per mode. For short and wide pumps (W ^> L), the sum 
over modes can be approximated by an integral. 

Similar results can be derived for the normal pump, 
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[4k%K 2 cos 2 (2k x L) + {k 2 x + Klf sin 2 (2A; x L)] 

Both the heavily doped contact (Qq~°°) and undoped 
contact (Qg~°) normal cases can be derived from the 
above by taking the x momentum in the leads K x to 
infinity or k x respectively. 



In Fig. [2] we represent the pumping response \q 01 
Eq. (§ as a function of q and kp in the two graphene 
setups, and compare these with the results of the normal 
pumps. For graphene we only show results for electron 
carriers. If the leads are heavily doped the pumped cur- 
rent reverts sign when the carriers in the central pump 
region are changed from electrons to holes [cf. Eq. 
In the case of undoped leads, the pumped current re- 
mains the same and only reverts sign if one also reverts 
the pumping cycle. 

In all four panels, the dashed line delineates the thresh- 
old between evanescent and propagating modes. The 
most evident feature in the graphene pump with heav- 
ily doped leads [Fig. [5|a)] is the contribution of evanes- 
cent modes to pumping close to the Dirac point. This 
contribution is absent in the other three cases, in which 
only propagating electrons are pumped. The reason for 
this is the high transmission of evanescent modes in 
graphene, which can be attributed to chirality and the 
Klein paradox 13 . Unlike in the normal pump with heavily 
doped leads, chirality conservation at the contact enables 
evanescent electrons to populate the graphene pumping 
region for modes within a window of width Aq ~ 1/L 
around q = 0. These evanescent modes contribute to 
pumping because they are sensitive to the onsite poten- 
tials Ui and have a finite amplitude at both contacts, 
so that charge transfer between them is possible over a 
pumping cycle. The finite contribution of the evanescent 
modes in graphene is in striking contrast to the vanishing 
contribution of electrons in the propagating mode q = 
(normal incidence), for which the transmission T g =o = 1 
is perfect at all energies because of the Klein paradox — 
these modes are therefore insensitive to driving and can- 
not be pumped. 

In a normal pump with heavily doped contacts [Fig. 
[2jb)], the large Fermi velocity mismatch suppresses the 
transparency of the contacts for all modes except those 
at resonant tunneling. This confinement creates sharply 
defined energy levels in the pump which are the origin of 
narrow regions of finite pumping, and results in a thresh- 
old kpL 1 below which no pumping occurs. In par- 
ticular, the contribution of evanescent modes to pump- 
ing is strongly suppressed (no contribution outside of the 
dashed line). Moreover, normal pumping in the limit of 
large lead doping is directed, meaning that for a given 
orientation of the driving cycle, the pumped current has 
the same sign for all energies. 

For graphene and normal pumps with undoped leads 
[Fig. [2|c) and (d)], all incoming modes remain propagat- 
ing in the pumping region. The main difference between 
graphene and normal pumps in this limit is the effect of 
the Klein paradox in the former, which suppresses pump- 
ing at q = 0, just as in the case of graphene with doped 
leads. Both pumps are open, and consequently there is no 
energy threshold for pumping. The sign of the pumped 
current is energy dependent, which is a generic feature 
of open pumps (including graphene pump with doped 
leads). 
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FIG. 3. (Color online) Total dimensionless response per 
mode x u f° r short and wide pumps (W ^> L), as a func- 
tion of the carrier concentration (parameterized by the Fermi 
momentum &f). Left panel: graphene pumps. Right panel: 
Normal-conducting pumps. Solid lines correspond to pumps 
with highly doped leads (which tend to 1/4 at /cpL ^> 1), 
while dashed lines correspond to pumps with undoped leads 
(which tend to 1/2). At the charge neutrality {UfL — > 0, see 
inset), evanescent modes allow for a finite charge transfer in 
the highly-doped-leads graphene pump, whose response satu- 
rates at a universal (sample-independent) value of 0.0288. 



In Fig. [3] we plot the total dimensionless pumping re- 
sponse X u (^f) of Eq. (|6| for the four types of pumps 
considered (left panel: graphene pumps; right panel: 
normal-conducting pumps). The regime of evanescent 
electron pumping in graphene with highly doped leads 
is visible as a saturation at kp = (see the solid curve 
in the inset). For wide and short pumps (W ^> L), this 
saturation value is sample independent and takes the di- 
mensionless value 
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0.0288, (8) 



which is the analogue of the minimal conductivity in the 
context of pumping } 1 5 * 16 i In contrast, all other pumps 
have a vanishing pumping response at depletion. At ener- 
gies kp > l/£, the pumping response rises to 1/2 and 1/4 
in the cases of undoped and highly doped leads, respec- 
tively. The normal pump with highly doped leads, how- 
ever, only operates above a finite carrier-concentration 
threshold corresponding to the position of the first reso- 
nant tunneling subband as mentioned above. 



IV. CONCLUDING REMARKS 



ing, to the extent that they provide the dominant contri- 
bution to the pumped charge when the system is operated 
close to the charge- neutrality point. Our comparison to 
results for normal pumps reveals that this effect is inti- 
mately related to chirality and the Klein paradox, and 
therefore arises from the unique low-energy properties of 
charge carriers in graphene. For the case of short and 
wide pumps, the evanescent pumping regime is charac- 
terized by a sample-independent universal value of the 
dimensionless pumping response. 

Practical considerations point towards additional ad- 
vantages of graphene-based quantum pumps. Firstly, in 
a realistic experimental setup, the principal pump driv- 
ing parameters are not the onsite energies Ui, but gate 
voltages Vi (see Fig. [I]) which manipulate the locally 
induced charge density rii(t) = Vi(t)Ci under each gate 
depending on the capacitances Q. The onsite energies 
then follow from the difference of the Fermi energy and 
the local position of the Dirac point. Neglecting details 
of the screening, for graphene ni = Wu? /(ttL), whilst for 
the normal case rii = Wuij (2i\V) (in the case of graphene 
with highly doped leads, a precise modeling would also 
have to take care of charge carriers populating the evanes- 
cent modes). In terms of these charge densities, the rel- 
evant response function is 



X n = X n c!et 



duj 
drin 



(9) 



involving the Jacobian between the u and the n vari- 
ables. Due to the divergence dui/drii oc n i close to 
the Dirac point, the experimental pumping response \ n 
for graphene pumps is expected to rapidly increase as one 
approaches charge neutrality, while it vanishes for normal 
pumps. Secondly, graphene-based quantum pumping 
promises to display an enhanced robustness against ther- 
mal effects. Thermal smearing of the pumping response^ 
occurs at temperatures of order £?l//cb, which are con- 
siderably higher in graphene than in normal pumps (this 
is also one of the reasons for the temperature robustness 
of other transport effects in graphene^). For the same 
reason, mechanisms limiting the pumping frequency are 
expected to be less stringent in graphene than in normal 
pumps. These considerations, together with the long co- 
herence length and high mobility in graphene, lead us to 
believe that graphene-based quantum pumps have good 
chances to become an experimental reality. 
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In summary, we find that evanescent modes can con- 
tribute significantly to graphene-based quantum pump- 
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